Macroscopic equations for the adiabatic piston 
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A simplified version of a classical problem in thermodynamics — the adiabatic piston — is dis- 
cussed in the framework of kinetic theory. We consider the limit of gases whose relaxation time 
is extremely fast so that the gases contained on the left and right chambers of the piston are al- 
ways in equilibrium (that is the molecules are uniformly distributed and their velocities obey the 
Maxwell-Boltzmann distribution) after any collision with the piston. Then by using kinetic theory 
we derive the collision statistics from which we obtain a set of ordinary differential equations for 
the evolution of the macroscopic observables (namely the piston average velocity and position, the 
velocity variance and the temperatures of the two compartments) . The dynamics of these equations 
is compared with simulations of an ideal gas and a microscopic model of gas settled to verify the 
assumptions used in the derivation. We show that the equations predict an evolution for the macro- 
scopic variables which catches the basic features of the problem. The results here presented recover 
those derived, using a different approach, by Gruber, Pache and Lesne in J. Stat. Phys. 108, 669 
(2002) and 112, 1177 (2003). 
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I. INTRODUCTION 

The so-called adiabatic piston is a long known prob- 
lem in classical thermodynamics, which can be stated as 
follows [l|, 0, y]. An isolated cylinder of length L, con- 
taining a gas, is divided by an adiabatic wall (no internal 
degrees of freedom), the piston, into two compartments 
(Fig. [T|). The initial condition is prepared in the fol- 
lowing way: the piston is kept fixed by a clamp at a 
given position Xq] the gases in the left (l) and right (r) 
compartments are in equilibrium defined by their pres- 
sure, temperature and volume: Pi,r,Ti,r,^i,r- By as- 
suming that the two gases are perfect and composed by 
Ni = Nr = N molecules with equal masses m, the gas 
state equation Pi^r^i^r = NTi^r holds in botli chambers 
(wliere the Boltzmann constant is set to unity by rescal- 
ing the temperatures). Being the piston adiabatic, the 
two subsystems are in equilibrium even if T/ =/= T^.. At 
t = the clamp is removed and the piston is free to move 
without friction with the cylinder. The nontrivial ques- 
tion is to predict the system evolution and the final posi- 
tion of the piston and of the thermodynamic quantities. 

At the beginning of last century, the above setup was 
used as experimental device for measuring the ratio Cp/cv 
of the specific heat of gases [Jl , that is linked to the pe- 
riod of the piston oscillations. Renewed interest on the 
problem has lead to recent experiments [a, m, ■ 

Meanwhile, several attempts were made to predict the 
final equilibrium state by using the laws of thermody- 
namics only, ending in controversial answers. A naive 



application of the first two laws of thermodynamics lead 
to the (wrong) conclusion that the equilibrium conditions 
is Pi/Ti = Pr/Tr-. A more careful treatment [3| shows 
that the correct answer is Pi = Pr- However, such a con- 
dition says nothing about the final position of the pis- 
ton and gas temperatures, which remain undetermined. 
Therefore, equilibrium thermodynamics cannot predict 
the final state. To shed light on the problem one has to 
cope with the non-equilibrium process that occurs after 
the clamp removal. 

From a microscopic point of view the adiabatic piston 
problem for ideal gases (non-interacting particles) can 
be described in terms of a one dimensional model, where 
the piston is a heavy particle of mass M much larger that 
the mass m of the gas molecules, which collide elastically 
with the piston. As argued by Feynman, the system first 
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FIG. 1: Sketch of the adiabatic piston. The subscripts I and 
r indicate the left and right compartments, respectively. 



converges toward a state of mechanical equilibrium with 
Pj. « P;, consistently with the thermodynamic predic- 
tion. Then, the pressure fluctuations, which are asym- 
metric because Ti ^ Tr, very slowly drive the system 
toward thermal equilibrium T^ = Ti [3|. In this way the 
final position of the piston and the thermodynamic quan- 
tities are determined. 

More recently, the problem was the subject of renewed 
attention, mainly stimulated by the talk of Lieb at the 
Statistical Physics Conference in 1998 Q, and by the 
connection of this problem with the ph ysic s of mcsoscopic 
systems [3, [lO] and Brownian motors [11| . 

Among the first attempts to understand quantitatively 
the time evolution of the adiabatic piston wc mention 
Crosignani et al. [l^l who introduced a set of ordinary dif- 
ferential equations for the macroscopic observables. How- 
ever, this model was only able to account for the position 
of the piston in the state of mechanical equilibrium and 
not for the final thermodynamic one. 

Remaining in the framework of ideal gases, a system- 
atic investigation in statistical mechanics terms, together 
with numerical simulations, has been carried on in the 
last decade by Gruber and coworkers [Ij, [Ij, [H, [I^, [l3, 
llSl . [ig . |2(| . In these works the problem has been ex- 
amined in several limits (see Ref. [3| for a review). In 
particular, in the thermodynamic limit taken by letting 
the system size L and the piston mass M to go to infinity 
by holding fixed the ratios po = N/L and R = mN/M, 
it has been shown that the system evolution can be re- 
duced to a set of ordinary differential equations for the 
macroscopic observables (i.e. the gas left/right temper- 
atures, and the moments of the piston velocity). These 
equations were obtained by using the Liouville and Boltz- 
mann equations. Within such an approach it is possible 
to control the deviations from Maxwell-Boltzmann distri- 
bution for the gas velocities, observed in the simulations, 
and a whole hierarchy of equations can be written for 
all moments of the piston velocity. Remarkably, these 
equations describe not only the reaching of mechanical 
equilibrium, which comes from the treatment at zero or- 
der in ^ [131, but also the final equilibrium state, which 
comes from the first order terms in ^ [1^ . These analyt- 
ical results have been shown by the same authors to be 
in agreement with numerical simulations of the ideal gas 
piston problem. Though it remains open the problem of 
a detailed description of the early stage of the dynamics 
in which the presence of shock waves has an important 
impact on the dynamics. Some recent attempts in this 
direction can be found in Ref. |2ll |. 

When the initial pressures are different, the system 
phenomenology can be described as follows [3|. In a 
first stage, the piston oscillates driven by the pressure 
difference. These oscillations are then damped till the 
"mechanical equilibrium" state, Pr ~ Pi but Tr ^ Ti, 
is reached. Then, as argued by Feynman, it follows a 
regime controlled by the asymmetry in the fluctuations 
felt by the left/right walls. This phase is characterized by 
a very slow approach to the thermodynamic equilibrium. 



Pr ~ Pi and Tr ~ Ti, with the piston position fluctuat- 
ing around the middle of the cylinder. In the oscillatory 
phase, both experiments |5|, [al , numerical computations 
and analytical arguments |17l.ll8l.l20l| have shown the exis- 
tence of two different regimes: weak and strong damping, 
the relevant parameter being R. For R < Re the adia- 
batic oscillations of the piston are weakly damped, while 
for R > Re they arc over-damped. Re being 0(1) [18j . 

Still in the context of ideal gases, it is worth mention- 
ing some recent approaches based on dynamical systems 
theory that have been developed by Chernov, Lebowitz 
and Sinai [23. In this context, also the case of gases 
starting from non-equilibrium conditions has been con- 
sidered [23, [2J] . Clearly, the ultimate goal would be to 
quantitatively understand the behavior of the system for 
an interacting gas, but this seems to be still too ambi- 
tious. Indeed only very few studies analyzed the case of 
gas composed by interacting particles [25|. 

In this paper, we consider a limiting case which has 
the advantage of being more tractable while displaying 
most of the non-trivial features of the problem. The ba- 
sic idea of our approach is to assume that the gases in the 
two compartments are composed of interacting molecules 
and thus characterized by a relaxation time toward the 
equilibrium state. Our main hypothesis is that this time 
is very short compared with all the other characteristic 
times of the system. In particular, we require that any 
fluctuation away from equilibrium (which is character- 
ized by homogeneously distributed gas molecules with a 
MaxwcU-Boltzmann velocity statistics) induced by the 
collision with the piston is re-adsorbed before the new 
collision with the piston walls. Physically speaking, the 
efficient re-adsorption of the fluctuations means that the 
(mechanical) work done by the piston is immediately con- 
verted into heat; an obvious consequence is that shock 
waves are ruled out. For the sake of simplicity, we also 
assume that the gases follow the perfect gas law. These 
hypothesis make the problem tractable while retaining 
the basic phenomenology of the original problem. 

Although a microscopic model of a gas able to fulfill 
the above requirements may sound rather artificial, at a 
practical level such a "microscopic model" can be easily 
implemented on a computer. The basic idea is to start 
with an equilibrium configuration with temperatures Ti^r 
for the gases, and then reinitialize the gas molecules as 
soon as one particle collides with the piston. The tem- 
peratures are recomputed after the collision and used for 
extracting a new configuration of the gas molecules. The 
procedure is then repeated. In the following we shall call 
such a model randomized gas. Even though, no actual in- 
teraction among the particles is actually considered, one 
can think that the re-generation of the gas configuration 
from an equilibrium one (but with the new temperature) 
is the result of such "unresolved" interactions. 

With the above assumptions for the gas, we will derive 
a set of ordinary differential equations for the time evo- 
lution of the macroscopic quantities describing the state 
of the system. Indeed the fact that the gas is always ho- 



iTLOgcncous and following the MaxwcU-Boltzmann distri- 
bution allows us to compute the joint probability density 
function that, in a given state of the system, the first 
colliding gas particle hits the piston in a time t and with 
a velocity v. Then, by averaging over this joint distri- 
bution the energy and momentum exchange due to the 
collisions with the piston, we derive the evolution of the 
macroscopic observables. The minimal set of variables 
required to have a closed set of equations is made up of 
the gas temperatures, the mean piston position, the first 
and second moment of the piston velocity. The second 
moment is required for accounting the piston fluctuations 
which, as argued by Feynman, are crucial for recovering 
the correct thermodynamic equilibrium [2, \lM, UmI • The 
equations are derived perturbatively up to the first order 
in ^. As we will see, though with a different approach 
and assumption, these equations are very similar to those 
derived by Gruber and coworkers [13, [l3 , ii^ particular 
at the zeroth order in ^ they are identical. 

We compare then the evolution of the system obtained 
by simulations of the ideal and randomized gas. In par- 
ticular, the agreement in the first (mechanical) regime is 
quantitatively perfect in the case of the randomized gas. 
While in the second regime, which is dominated by the 
fluctuations, the agreement seems to be only qualitative. 
Somehow surprisingly, we found that, in this regime, a 
better quantitative agreement seems to be possible dis- 
regarding some 0{jj) terms. However, with such terms 
excluded the equipartition of energy at equilibrium is vi- 
olated by the piston. Some hints to explain these findings 
could come by higher orders terms in the expansion. Un- 
fortunately, the computation of the higher order terms is 
very cumbersome. Since the newest aspects of our work 
is in the proposed derivation and in the introduction of 
the randomized gas model, we present in this paper the 
all approach up to first order in ■^. 

The paper is organized as follows. In Sect. |TT] we 
present our approach based on the collision statistics and 
derive the equations for the macroscopic observables. In 
Sect, mil we compare the results of the model with those 
obtained by simulations. Discussions and conclusions can 
be found in Section llVl In order to avoid long appendices, 
the technical material, with the detailed derivation and 
all the formulas needed to make explicit the equations, is 
presented as electronic supplementary material (26[. 



II. DERIVATION OF THE MACROSCOPIC 
EQUATIONS 

The underlying idea of our approach is to derive a set of 
deterministic dynamical equations for the "macroscopic" 
variables describing the evolution of the thermodynamic 
state of the system under the assumption that, at any 
time, the gases in both chambers are perfect and at equi- 
librium. In other words the gases are able to instanta- 
neously dissipate the fluctuations induced by the colli- 
sions with the piston. Thus a MaxwcU-Boltzmann equi- 



librium state holds always in both compartments but, in 
general, with different temperatures and volumes. 

While the above hypothesis define the macroscopic 
state of the gas, for the piston the problem is more sub- 
tle. We would like to describe its motion on times longer 
than the single collisions, that is to average its instan- 
taneous position and velocity {X, V) over the collisions 
so to obtain a deterministic (macroscopic) trajectory de- 
fined by the average position x = X and velocity Vx = V 
The symbol [. . . ] denotes the average over the collisions. 
As discussed in the introduction, it is crucial to account 
also for the fluctuations of the piston velocity. For this 
reason, the second moment of the piston velocity V^ is 
included in the description. 

In the thermodynamic limit we will consider, one can 
argue that the fluctuation of the piston position can be 
safely ignored. This means that in the following we will 
consider the piston position as a deterministic quantity 
and we shall use only the mean piston position x. 

Given the piston position, the gas is characterized by 
the temperature Ti^r-, and volume fli^r, with 0,1 = x and 
ilr = L — X (we assume a Id geometry for the sake of 
simplicity). Being perfect gases, the pressures are given 
by the equation of state Pi^r^i^r = NTi^r- 

In the sequel we show how to derive a_set of differential 
equations for the evolution of x, v^, V"^ and Ti^r- In 
order to keep the presentation as simple as possible, here 
we shall sketch how the equations can be derived and 
the averages performed skipping all the algebra of the 
computation, which is detailed in [26|. 



A. Macroscopic equations from the collision rule 

For the formal derivation of the deterministic equa- 
tions, we only need the above discussed assumptions and 
a microscopic ingredient: the elastic collision rules 
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Primes denote postcollisional velocities, and v the collid- 
ing gas particle velocity. The quantities we are interested 
in are the time derivatives of the macroscopic observables 
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where we set the Boltzmann constant fc^ = 1. The time 
derivatives should be computed starting from the colli- 
sion rules as the averages ([. . .]);,r = [•■•]; rl^^ suggest. 



St being the mean collisions time (for a more precise and 
operative definition see Sect III Cp . The subscripts I, r de- 
note averages performed over the collisions with particles 
residing on the left/right compartments. 

V'^ — V which evolves 
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where we used Eqs. ([3|) and fl]). It should be noted that, 
at this level, the piston is completely described by Vx and 
(jy . This amounts to assume that its velocity distribution 
is Gaussian 
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Since at the initial time t = 0, one starts with Vx ~ Q and 
Gy = 0, the above probability distribution is initially a 5 
function. Plugging ([T]) into ((3]l5|), we obtain 
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The equation for <Ty ^ is obtained from ([6]) and ([4]) by 
using the collision rules ([T]) . In ([8|) and ^ , the prefactor 
N appears as a result of a time rescaling, that sets the 
time unit to the average collision time, which is order 
1/A^. Said differently, the change of the gas temperatures 
due to the collision with the piston is order 1/A^. 

Notice that, for reasons that will become clear in the 
following, the average of the type {V v) is different from 
Vx{v). We anticipate that this difference is not due to 
a breakdown of molecular chaos hypothesis (as one may 



naively think) but to the fact that the collision statistics 
depends on the instantaneous value of the piston velocity. 
More explicitly, Vx{v) represents only the zeroth order 
term of (Vv), and a term coming from the fact that V 
is a fluctuating quantity will also appear. 

Notice also that the above equations conserve the total 
(gas plus piston) energy 
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The consistent (first order in -^ ) equations can then be 
obtained from ((Sl llip by expanding the various prefactors, 
performing the limit N —> oo and suitably expanding 
around it. Since the procedure is delicate, we proceed 
step by step. 

B. Thermodynamic limit and formal expansion of 
the equations 

First of all we have to specify the limiting procedure, 
that as explained in [3| can be done in different ways. 
We are interested in the limit N, A/, i ^ oo in which 
we keep fixed po = N/L and the nondimensional mass 
ratio R = Nm/M . We have now to expand around this 
limit retaining all terms which are first order in ^ (and 
consequently at the first order in l/N). Aiming to make 
explicit the zeroth and first order terms we formally write 
the averages as: 
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where the first and second terms on the r.h.s are the 
zeroth and first order terms of the expansion. How to 
explicitly perform such an expansion will be explained in 
the following subsections. We warn the reader that for 
maintaining the notation as compact and explicit as pos- 
sible in the following we adopt the convection to indicate 
with ([. . .J)'-^^ all averages which are 0{^) irregardless 
if this comes from the expansion of the average or from 
the averaged quantity. For instance, by direct inspection 
of Eq. ^ at equilibrium, one easily realizes that ay is 
0{^). Therefore, we shall always indicate its average 
with (cTy)'^^. Finally, notice also that all the terms in- 
volving powers of {V — Vx) vanish at the zeroth order. 
Keeping in mind these simplifications, the (expanded) 
equations become: 
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Before sketching the way the above averages can be com- 
puted (see Sect. Ill CI and [2^) we briefly discuss some 
properties of the above equations. 



The first observation is that Eqs. (|14m7p ensure the en- 
ergy conservation p^ both at the zeroth and first order, 
meaning that the expansion is consistent. Notice also 
that the relative importance of the various terms is not 
the same at all times. As the system evolves, their rela- 
tive weights change corresponding to the different stages 
of the evolution, briefly summarized in the introduction 
and detailed in the following. For example, consider ay. 
At the beginning ay — 0, then it grows until it reaches 
its equilibrium value. Consequently, the terms which in- 
volve the velocity fluctuations are not important at the 
beginning. While they become 0{^) and play a crucial 
role in the final stage of the system evolution. The op- 
posite is true for the terms involving the average drift 
Vx, which is close to zero in the second (Brownian) stage 
of the evolution, and large in the first (mechanical) part 
of the system evolution. Among the terms involving the 
piston velocity fluctuations, we should mention a special 
role played by those in which it does not (explicitly) ap- 
pear ay. These are the terms of the form {{V — Vx)vx) l ^ 
and {{V~Vx)v)l I . as we will see they will be both propor- 
tional to ay and, as anticipated, find their origin in the 
way the fiuctuations of V affect the collision statistics (see 
next subsections and the supplementary material [26|). 

However, a closer inspection shows that {{V — Vx)vx)il 
is very small at all times. In the first stage the fluctua- 
tions are negligible while in the second stage the average 
drift is very small. Though we retained this term in the 
equations and in the numerical simulations, one can show 
that they can be removed without problem. Differently 
the terms {{V — Vx)v)i^ are very important in the fi- 
nal stage of the evolution and, as discussed below, for 
obtaining the correct equilibrium state. 



We are still left with performing the infinite volume 
limit. We anticipate here that all the terms appearing in 
the averages arc proportional to either l/x, when coming 
from a left-chamber average, or 1/(L — x) when coming 
from a right-chamber average. There is no other depen- 
dence on a; or L in equations ([2115]) and, consequently, 
Eqs. (|14II17[) . This implies that one can simply intro- 
duce the hydrodynamic time tn = tjL and the rescaled 
X coordinate xu ~ xjL. With this rescaling, L does not 
appear anymore in the equations and there is no need to 
perform the limit, meaning that there are no corrections 
to the results due to finite size effects. For a simpler 
comparison with the simulations, we will keep writing 
in the following the equations for finite values of L; the 
corresponding expressions in the hydrodynamic time can 
be simply obtained with the above substitution, that in 
practice corresponds to set L = 1. 



1. Equation at the 0*'' order in jj: Mechanical Regime 



Let us start a closer inspection of the equations start- 
ing from the zeroth order terms, i.e. assuming -^ ^^ 0. 
The velocity fiuctuations of the piston are ignored (mean- 
ing the motion of the piston is purely deterministic) and 
the final equilibrium position depends on the initial con- 
ditions. The result is anyway nontrivial. As discussed 
in [i3, ll3 J having considered the dynamics and not only 
the thermostatics (which only tells us the equality of pres- 
sures), we can now determine the mechanical equilibrium 
position. In the sequel, we shall show that at the zeroth 
order in jj we obtain, by using a different approach, the 
same equations of Gruber and coworkers [l8|, and fol- 
lowing them we sketch how the mechanical equilibrium 
point can be computed. 

For obtaining the equations at the 0*'* order, we need 
to set -Ij = and to ignore all averages indicated with 
the superscript ^^^ in (|14p . ([TBI) ^^^^ (|17p . In other words 
wc only need the averages 
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See the supplements [26| for the derivation of the above 
expressions. Assuming Vx ^ y/Ti^r/'m (which is reason- 
able for realistic values of the physical parameters), and 
expanding ([T51) and ([TO]) in Vx, the equations ([Til) ^-^^^ 
(ITBlHTl) reads 
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where for the friction coefficients it holds 7 = 7; + 7r 
ensuring energy conservation p2p at the 0*'' order; from 
(|18m9p . at the lowest order in Vx one has: 
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we remind that Qi = x and fi,. = L — x. Note that 
apparently, in the limit L —^ 00, 7; ,. -^ 0, this is not 
the case if the hydrodynamic rescaling is properly ap- 
plied. Indeed, taking the hydrodynamic limit the above 



expression remains unchanged, keeping in mind that in 
this case the " hydrodynamic volumes" are flf = xu and 
il^ = 1 — xh- In particular, the damping coefficients 
go to a finite value also in the hydrodynamic limit. It 
is worth remarking, that the linearized equations (|20ll22p 
coincide with those derived in Ref. ^M, with a different 
method. In the absence of friction, one can easily see 
that they describe a purely adiabatic transformation of a 
one-dimensional perfect (mono-atomic) gas. Indeed the 
first term on the r.h.s. of Eqs. (pn)) is simply the pressure 
difference on the two sides of the piston, while the first 
term of Eq. (PT|) and (P^ can be obtained by differenti- 
ating with respect to time the equation of an iso-entropic 
process, namely 






(24) 



where Cp/c^ = 3 is the specific heats ratio and the initial 

conditions fix C(,r = T/ ,.(0)ri^'J, (0). In the absence 

of the friction terms this would give rise to periodic oscil- 
lations of the piston. As discussed in [13 . [l8j , the friction 
terms are responsible for the irreversible evolution toward 
a state of mechanical equilibrium for which 
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where the tilde indicates the mechanical equilibrium 
quantities. Notice that in this framework irreversibil- 
ity naturally emerges as a result of the averaging over 
the collisions [231. Eq. ((25)) is also the result of ther- 
mostatics, but it is not enough to determine x and T; ,,■ 
Indeed; the full dynamics given by (|20ll22p is needed to 
predict such mechanical equilibrium point, as shown in 
the sequel where we brieflysummarizc the results first 
derived in Gruber et al. [i8[. First, notice that ((25|) 
together with ^ tell us that fij = 2Toni,r/L, with 
To = E/N = {Ti{0) +TriO))/2. So that' the equi- 
librium pressure is P == 2NTo/L. Second, defining 
Z = ^/Tix— ^/T\-{L — x) and by using (|20ll22p one can eas- 
ily see that dZ/dt = 0{v^), which if Wj, <C 1 means that 
Z is conserved. Z = const provides the missing condition 
to determine the mechanical equilibrium. The resulting 
equation for the equilibrium point is therefore |18j : 



Vx)vx) ) allows for energy exchange among the two com- 
partments mediated by the fluctuation of the piston. Of 
course, such terms start to play a role once the mechani- 
cal regime (described by the 0*'' order terms) is finished, 
i.e. when the fluctuations of the piston become relevant. 
This regime driven by the fluctuations results from the 
expansions in -^ and 1/A'^ which, as it happens com- 
monly in Brownian motor-like systems |ll| . are inter- 
twined and add new (sometimes unexpected) features to 
the dynamics. In particular, in this case one can show 
that Eq. (|14m7p evolve toward a nontrivial stable fixed 
point corresponding to the thermodynamic equilibrium, 
i.e.: 
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Ti = T, = T, 
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eq 
M 



(27) 



The last equality cannot be explicitly seen from ([1| 
which simply states that at equilibrium (fx^)'^^) = 
™(„2^(o) _,_ ^(^ „ y^)y)ii). As we discussed, {{V - 

Vx)v)^^' oc ay and with the explicit computation at equi- 
librium [2y] one can see that cry = Teq/M, which is a 
pleasant result since it is in agreement with the condition 
of equipartition of energy. Notice that Eq. (P7)) suggests 
to interpret May ~ Tp as the temperature of the piston. 

We conclude this subsection mentioning that the above 
equations are similar with the ones obtained by Gruber, 
Pache and Lesne [19[. Due to the very long expressions 
involved in the equations at the first order, we could not 
decipher whether they are exactly equal. At the end of 
next section we shall discuss the possible source of dif- 
ferences. However, we stress that in these two works a 
different approach and different assumptions were made 
on the gases. In particular, Gruber, Pache and Lesne 
derived the equations from an expansion of the Boltz- 
mann and LiouvuUe equations. In some sense, our and 
their different assumptions can be seen as two different 
ways to close the hierarchy of equations to the second 
order and one should expect the phenomenology of the 
two equations to be, at least, qualitatively similar. 

Section IIIII is devoted to compare the evolution of the 
macroscopic observables obtained by integrating p4j|17p 
with that of the microscopic model. In the following sub- 
section we detail the procedure by which the averages can 
be computed. 



'TiS; - ^jTr{L -i)^ y/Td^x{Q) - yjTr{Q){L - x{0)) , 

(26) 
which should be solved for x after plugging T; = 2T{)x/ L 
andfr =2To(l-i/L). 

2. Equation at the 1^^ order in ^: Brownian Regime 

When the first order (in -jj) terms are retained, the 
terms in ay (among which, as discussed above, we have 
also to consider the terms {{V — Vx)vY^^ and {{V — 



C. Explicitation of the averages 

In order to finalize our program we have now to make 
explicit the averages in (|14m7p . Let us start by mak- 
ing explicit the formal expression of averages such as 
([. . . ])/,,, which should be interpreted as follows. Denote 
with Gi^r{t,v\V) the probability of having a left/right 
collision in a time t with a velocity v for the r/1-particle 
conditioned to a realization in which the piston has veloc- 
ity 1/, and indicate with G their sum, i.e. G = Gi + Gr, 
which is normalized to 1. Then performing the average 



{A{v, V, t))i^r of a generic funetion A{v, V, t) means: 

/ dt JJ dvdVP{V)Gi^r{t, v\V) A{v, V, t) 



{A{v,V,t))l,r^ 



J dt JJ dwdy t P[V)G{t, v\V) 



(28) 

where St = J dt JJ dwdV^ t P{V)G{t, v\V) is the mean col- 
hsion thne. 



1. Derivation of Gi^r{t,v\V) 

Let us now derive Gi and Gr- In particular, we shaU 
compute them in the thermodynamic hmit by exphcitly 
considering the terms order 1/A^ and consequently ^. 
These terms are those entering the averages we indicated 
with ([...])«. 

We start by the equilibrium distribution of the gases, 
which is uniform in the particle positions y and Maxwell- 
Boltzmann for the velocities v. 



Pi{y,v) = Pi{y)Pi{v) ^ .- 6 -^1 



Pr{y,v) = Pr{v)Pr{v) = 



x^/2wrrFWi 
1 



(29) 



[L - x) ^2Tnn-^Tr 



that, under our assumptions, describe the gases at all 
times. Note that the above distributions depend (para- 
metrically) on the dynamical variables x and T/ ,,. 

From the equilibrium joint distributions p9p we can 
derive the probability density gi^r(t\v,V) that a particle 
on the left/right collides in a time t given its velocity 
V and the macroscopic state of the system defined by 
the temperatures, the piston position x and velocity V 
and the equation of state. Of course, most of the weight 
to such an hitting probability comes from particles that 
are close to the piston and that have a large (negative) 
relative velocity with it. These particles are far from the 
bulk of the gas and in this derivation we assume that we 
can use the evolution of free particles to compute their 
hitting time with the piston. Within such an assumption, 
it is easy to realize that the probability gi.r{y, v) is simply 
obtained through a change of variables from (P5)) 



9iit\v,V) 



e{t)e ( ^TT-t) e{v~v)^^-^ (30) 



v-V 



e{t)Q{v-v) 



v-V 



g^{t\v,v) = eit)e(^-t)eiv-v)y-^ (3i) 

\V—v J L—x 



Q{t)Q{V~v) 



V-v 
L — x 



<d being the unitary step function. Note that in the sec- 
ond expression we ignored the second Q function. This 
is justified by the fact that t indicates the time be- 
tween two consecutive collisions, which for iV ^ 1 is 
always much shorter than the time needed for a parti- 
cle to travel along a whole chamber. Further, consider- 
ing only positive times, the functions gi^r{tiV\V) do not 



depend on t. We shall then use the compact notation 
fi.r{v\V) = gi_r{ti'v\V). Notice also that we use the in- 
stantaneous piston velocity V and not its average v^, 
while we ignore the fluctuations in position. Moreover, 
we neglected possible correlations between the velocities 
and positions, which amounts to implicitly assume molec- 
ular chaos. 

Given gi^r{t\v^V)^ the joint probability of having an 
impact of a particle with v in a time t is : 



5i,r(i, v\V) = gi,r{t\v, V)pi^r{v) 



(32) 



The probabilities Fi,r{tm) that a left/right particle col- 
lides in a time < t„i is then given by 

FiAtrn)= dt dv giAt,v\V)^trn dv fiAv\V) (33) 



(34) 
(35) 



which we rewrite as 

Pl,r{tm) — tmHl,r{x,V,TiA . 

with 



Hi{x,V,Ti) = 

Hr{x,V,Tr)'- 



-e ^^1 



dv- , 

V \/2'Km^^Tix 

V 

dv- 



V-v 



\/2Trm ^Tr{L — x) 



In the following we will use the shorthand notation Hir = 
Hi^rix,V,TiA and hi^r = Hi^x^Vx^Ti^. 

Considering that we have N particles on both the left 
and right, the probability densities that one of them on 
the left/right impacts the piston in x with a velocity v 
are given by 

Gi{t„,,v)^N[l-Fi{t,r,)f-\l~Fr{t,^)ffi{v\V) 



Gr{tra,v)^N[l-Fr{tm)r~'[l~Fl{tr,-.Tfr{v\V). (36) 

We can now perform the limit N —^oo and tm — > holding 
Ntm ~ T fixed. Noticing that Fi.r{tm) -^ THi,r/N and 
expanding [1-Fi^r]^ =exp[N\n{l-Fi^r)] « exp[7V(-F/.^- 
F^r/'^)] ^^'^ retaining terms only up to 1/A'', we find 



Gl,r{T,v\V)^ 

with 



J-H — TTT T 



N 



2N 



Gl,r{T.v\V) [37) 



GiAr,v\V) ^ e--^"^+"^^ fiAv\V) . 



(38) 



By recalling that N = AIR/m, and aiming to retain only 
the first order -^ terms (|37l) can be rewritten as 



GiAr,v\V) = GiAr,v\V) + 



m f hirT {hf + h'^) 2 
M V~R 2R ^ 



(39) 

GLriT,v\Vx) 



where we substituted h in place of H and v^ in place of 
V in the part which is already at the first order in -^j-. 
However another contribution to the 0{jj) term comes 
from the expansion of the first term in 



2. Expansion of the average 



We need now to further expand (|39p in V around Vx, 
this can be accomphshcd by Taylor expanding G 

G{t,v\V) = GiT,v\vx)+dvG{T,v\vx)iV^Vx) 

+ ld^G{T,vK){V-v,f. (40) 

Noticing that JdVP{V){V - v^) = and that 
J dVP{V){V — VxY = cTy ■^e can write, for example, 
the (expanded) average collision time: 



5t 



1 



hi + hr 



m hf + h^ 

Tl R{hl + hrf 



(41) 



hl±hL^l 1 1 dtdvtdlG{T,v\vx) 



Finally, we can write the correct expansion of ((28|) at 
the first order in ^, by plugging all the expanded terms 
in (P5|) to have the zero and first order terms of the aver- 
age of a generic observable, {A)r.i = (A)], / + {A)l j . The 
result is: 



(A) 



(0) _, 



(^)l!.^ 



(hi+hr) dTdvdVP{V)Gur(T,v\Vx)A (42) 

dTdvdVP{V) dvGi^TAvxW ~vx) 

m hir 2 



+ ^d^GiATMv.)(y-Vx)'' 



A 



dtdvtd^G{T,v\vx) 



M R{hi + hr) ^ 
dv//,.(w|«.)A(43) 



{hi + hr 



which can be used to compute all the averages in (|14t - 
I17p . However we mention that there are exceptions to the 
above recipe. For instance, the average of A = {V — VxY , 
whose result is (cr^)'-^-' and not (ct^)^°^. This is due to 
the convention adopted for the average: remember that 
we chose to write the superscript '^' also when the jj 
order comes from the averaged quantity and not from 
the expansion of the collision distribution, like in this 
case. 

Notice that in the first term of the r.h.s of (j43|l . as for 
the average collision time, we expanded G with the aid 
of ((40l) . Depending on the observable A, which may have 
or not a linear term in V — Vx, also the first derivative of 
G may appear. 



III. COMPARISON BETWEEN MODEL AND 
MOLECULAR DYNAMICS SIMULATIONS 

In this section we compare the evolution of the macro- 
scopic observables given by (|14m7p with numerical simu- 
lations of the microscopic model. We consider two kinds 
of microscopic simulations: the ideal and the random- 
ized gas. The latter is meant to fulfill the assumptions 



we made on the gas in deriving the equations. Let us 
now better clarify how this is realized. 

From a computational point of view, it is very easy to 
realize the randomized gas, the idea is to let the gases re- 
lax in an artificial manner through a randomization pro- 
cedure. More precisely, the simulations are performed in 
the following way: we generate an equilibrium configura- 
tion of the system, corresponding to given values of the 
macroscopic observables. Then the system is let to evolve 
up to the first collision of the piston without interactions 
among the gas particles. After the collision, the energy 
of the gas containing the colliding particle changes; the 
fast relaxation of the system is mimicked by updating 
the gas temperature (corresponding to the new energy) 
and redrawing an equilibrium configuration of the gases 
corresponding to the new temperatures (and volumes). 
Then the process is iterated. Notice that in this way one 
keeps track of all the observables except the piston ve- 
locity fluctuation, which needs some kind of average to 
be defined/measured. Collisions are evaluated, as for the 
ideal gas, with an event driven algorithm. 

In Figure [2] we show the evolution of the piston posi- 
tion by numerical simulations of the ideal and random- 
ized gas for two different values of R (namely i? = 10 
(top) and R = 0.2 (bottom). As discussed in the in- 
troduction (see [3, [13 for a more detailed treatment) 
these two choices correspond to the case of strongly and 
weakly damped oscillations, respectively. In both cases 
one has two regimes: mechanical and Brownian. The 
former is characterized by damped oscillations of the pis- 
ton, which in the strongly damped case are very few or 
nonexistent. For the ideal gas, as discussed by Gruber et 
al. 0, [13, 13 the detailed damping of the oscillations de- 
pends on the presence of shock waves, which are absent in 
the randomized gas. In fact, the latter is damped much 
more efficiently than the former. As shown in the bot- 
tom panel, for R = 0.2 the evolution is weakly damped 
and many oscillations are observable. In this case, the 
piston is very slow and the gases perform quasi-adiabatic 
oscillations, which are damped more and more mildly as 
i? ^ 0. Again, also in this case the damping appear to 
be much more efficient in the randomized case. 

The Brownian motor- like regime [ll| , occurs when the 
oscillations are completed damped, i.e. the mechanical 
equilibrium is realized with approximately equal pres- 
sures Pi ^ Pr, which differ only for terms 0{^). From 
now on, both the ideal and randomized gas remain in 
a state of marginal equilibrium approximately along the 
isobar Ti/x = Tr/{L — x), which is the prediction of ther- 
modynamics (shown in Fig.[3]for the ideal gas only). The 
Brownian motor-like mechanism is responsible for heat 
transfer from the warmer to the colder chamber medi- 
ated via the piston fluctuations [111, [2^ . This stage oc- 
curs on a long timescale (proportional to M). For re- 
alistic values of the various parameters, one can realize 
that the timescales necessary for reaching the thermo- 
dynamic equilibrium state are enormous. It is therefore 
very difficult to observe this regime in experiments, and 
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FIG. 2: (Color online) (top) Comparison between the evolution of the piston position x(t) in a simulation of an ideal gas 
(red solid line) and of the randomized model (blue, open circles). In both models we set A'^ = 1000, M — 100 and L — 2000 
corresponding to _R = 10, at time zero T; — 40, Tr = 60 and x = 0.61/. To have a cleaned curve we performed an average of 
about 100 independent realizations, (bottom) The of (top) with N = 1000, M — 5000 and L = 2000 corresponding to i? = 0.2, 
and the initial state set as Ti = 150, Tr = 50 and x — 0.6L. Here, for the ideal gas, the oscillations last for a much longer time 
and are damped very slowly before the Brownian motor like regime sets in; while the randomized gas is much more efficient in 
damping the oscillations. 



well controlled numerical simulations are mandatory. In 
the final state the two chambers have the same temper- 
atures and the piston fluctuates reaching equipartition 
with the same temperature. It is worth noticing that for 
the ideal gas, since the molecules interact only through 
the collisions with the piston, the probability distribu- 
tion of the velocities of the gas molecules may (and ac- 
tually do) deviate sensibly from the Maxwell-Boltzmann 
distribution that is recovered only after the reaching of 
the final equilibrium state [3|, ll9[ . By construction, this 
problem is not present in the randomized gas which is 
forced to remain Maxwell-Boltzmann. 

In the following wc shall compare the evolution of the 
macroscopic equations with the simulations. 



A. Mechanical regime 



X, Ti^r and Vx — 0. Linearizing (|21j|22p one can easily 
recognize that the equation define an iso-entropic pro- 
cess, i.e. (p4)) . meaning that Ti{x) = Ti{x/xY^^'^'-'~^ and 
Tr{x) = fr[{L-S:)/{L-x)Y'-'''--^, which plugged in ^ 
and expanding \a 5x = x — x lead to the equation of a 
damped oscillator: 



qV5x R A5x 
di^ m dt 



R 



PL 



mN Cy x{L — x) 



5x^Q, (44) 



where we recall that Cp/cy — 3 and P — 2NTo/L (see 
Scc IIIBT]) . Since the friction coefficient does not modify 
the period one immediately gets: 



r = 27r. 



INmx{L - x) 
R{cp/cy)PL 



(45) 



Wc start the comparison by considering the oscillatory 
regime in the weakly damped case. Here we know from 
previous studies [l^ that the model is able to correctly 
predict the period of the oscillations. 

The period of the oscillations in the model can be es- 
timated from the linearized dynamics at the zeroth or- 
der in ^, Eqs. pl^ . As discussed in Sect. HTbI the 
oscillations will occur around a mechanical equilibrium 
position X defined by Eq. ([2S1). Then to recover the 
period of the oscillations it is enough to linearize (|201 
around the mechanical equilibrium state defined by 



this formula is at the basis of the measurement of the 
specific heat ratio in experiments [J| . 

In Fig. m we compare the ideal gas simulations with 
that of the randomized gas and the numerical integration 
of the determinist equations (|14m7p . As one can see, the 
ideal and randomized gas have the same period, while the 
damping is different, and the model is in perfect agree- 
ment with the randomized gas simulations. Integrating 
the equation at the zeroth order, we find x ~ 6.39 x 10'^ 
which agrees with the predicted value (pS)) and Eq. ([^S]) 
predicts T = 2.17 x 10^ which is in very good agreement 
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FIG. 3; (Color online) Evolution of the system in the tem- 
perature volume (TQ.) plane. The dotted line indicates the 
adiabatic phase T oc V~^ , while the solid straight line the iso- 
bar that characterizes the final phase (Brownian motor like 
regime). The rightmost arrow indicates the initial state of 
the gas in left compartment Ti{0) = 150 ^((0) = a;(0) = 0.6L 
which evolves (red dots) toward the final equilibrium in the 
middle. The leftmost arrow indicate the initial state of the 
gas in left compartment ri(0) = 150 0^(0) = L - x(0) = 0.4L 
which evolves (blue dots) toward the final equilibrium in the 
middle. The simulation is done with the ideal gas, for which 
the initial oscillations are much more evident that for the ran- 
domized one. 



with the the measured period in both the randomized 
and ideal gas. 

We mention that in Ref. [l^| a detailed study of the pe- 
riod of the oscillations was reported. Since our linearized 
equations coincide with those of Ref. [ig , we shall not re- 
peat here this study. It is however interesting to compare 
the first stage of the evolution of the ideal and random- 
ized gas with the model in the case of strong damping. 
In Fig.[5]we show the leftmost plot of Fig. [5] (top), as one 
can see although the model is unable to reproduce the 
oscillation of the ideal gas, its evolution coincides with 
that of the randomized gas. 

As discussed in Ref. @] , estimating the decay rate of the 
oscillations for the ideal gas is a nontrivial task, since it 
requires a detailed study of the dissipation mechanisms 
of the shock waves created by the piston motion. These 
shock waves survive for a long time in the ideal gas, while 
they lifetime is expected to be shorter in the presence of 
interactions, which should be able to decrease their coher- 
ence. The randomized gas represents a sort of limiting 
case of interaction in which the shock waves are com- 
pletely absent. Most likely the absence of shock waves is 
at the origin of the faster damping of the oscillations in 
the randomized gas and of the very good agreement be- 
tween its evolution and the one obtained from the macro- 
scopic equations. It would be interesting to investigate 
the transition between weak and strong damping in the 
randomized model. This is far from the aim of the present 




0.55 



0.45 



5x10"^ lOxlO'* 15x10"* 20x1 0'^ 



70 




















' 






60 
50 


\k 


* 




« , ffl g, 


□ D 


□ n n J 
^— n — 


40 


v ^ 


I 


^ 


¥ °tf ^° w 




□ □ n □ - 


30 


























70 
60 














11 


b 


^ 






□ ° n □ - 

□ ° an 


50 
40 


/y 




♦ ' 


W ^ °rf *° 




30 



























5x1 0^* lOxlO'^ 15x10"^ 20x1 0^* 

t 



FIG. 4: (Color online) Comparison between the simulations 
of the non-interacting ID gas, the randomized one and the 
prediction given by the macroscopic equations. The parame- 
ters are: Ti = 40, Tr = 60, L = 12000, xo/L = 0.6, M = 10^, 
A^ — 10"^ corresponding to R = 0.01. (top) Evolution of 
the piston position, (bottom) Evolution of the temperatures. 
Blue open squares refer to the ideal gas and red filled circles 
to the randomized gas, the solid line is the model p4ll7p . 



paper, however the above results suggest that the criti- 
cal value of R for the transition will probably be smaller 
(but of the same order) of that of the ideal gas which is 
order unity. 



B. Brownian motor-like regime 

As the mechanical equilibrium is reached Pi ~ Pr, 
the system slowly evolves, along an approximate isobar, 
driven by fluctuations toward the thermodynamieal equi- 
librium defined by Eq. (P7|) . which is the (stable) fixed 
point of the ordinary differential equations (|14m7|) . 

We shall now compare the evolution given by the 
macroscopic equations with that of the ideal and ran- 
domized gases in the Brownian regime. To minimize 
the possible differences between the ideal and random- 
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FIG. 5: (Color online) The same as the leftmost plot of 
Fig. EJtop) with superimposed the thick black line superim- 
posing on randomized gas data (blue empty circles) is ob- 
tained by the numerical integration of (|14m7p . The agree- 
ment of the model with the randomized gas is perfect also for 
the temperatures (not shown). 
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FIG. 6: (Color online) Comparison between the molecular dy- 
namics of a non-interacting ID gas and the model in the Brow- 
nian motor-like regime. Simulations have been performed set- 
ting the initial state as T; = 40, Tr = 60 and a:(0) = 0.4L with 
M = 500 and L = 12N with A^ = 10*: red triangles refer to 
the randomized gas and the blue squares to the ideal gas, 
respectively. The dashed line is the prediction of the macro- 
scopic equation, the solid line (which perfectly superimposes 
on the randomized gas data) is explained in the text. The 
simulation data are obtained by performing an average over 
about 10 realizations to reduce the fluctuations. 



izcd gases we performed a simulation which starts in a 
mechanical equilibrium state having the gas molecules 
distributed according to the Maxwell-Boltzmann distri- 
bution. For the ideal gas, this might be not the typ- 
ical situation: usually when the system arrives to the 
mechanical equilibrium from a non-equilibrium state, it 
may be strongly non-Maxwellian [3, [l^ . As exemplified 



in Fig.[6l where we show the deviation from the final equi- 
librium of the piston position, \x — Xeq\/L, the relaxation 
is exponential. The randomized gas relaxes faster than 
the ideal one, likely because the latter develops slightly 
non-Maxwellian distribution (we observed that the dif- 
ference in the relaxation time tends to diminish as the 
number of particles is increased). Note that the simula- 
tion has been performed with R > 1 because the time 
scale for reaching the equilibrium is controlled by M; 
with N = 10^, as here, working with R < 1 would have 
implied the necessity to reach too large time scales to 
study the relaxation. 

As shown in the figure, the macroscopic equations 
(dashed line) predict a relaxation slower than both the 
randomized and ideal gases. This came out as a surprise 
for us, because we were encouraged by the very good 
agreement in the mechanical regime, discussed in the 
previous section. With the aim of understanding such 
a difference we examined all the terms appearing in the 
macroscopic equations and we realized that the mismatch 
in the relaxation was due to the terms {{V — Vx)v)i^ 
that, as discussed in Sect. IIIB2| are those which en- 
sures the equipartition of energy at equilibrium, i.e. the 



fact that Mai 

V eq 



Teq. In particular, eliminating such 
terms from (|14lll7p (note that this is not affecting the 
conservation of energy), we found a perfect agreement 
between the macroscopic equations, now modified, and 
the randomized gas, as shown by the solid line in the fig- 
ure. Nevertheless, with such a modification the energy 
equipartition is not verified anymore and at equilibrium 



M<yVeq 



2T, 



eq- 



One may think that the agreement is incidental, and 
so we did a more severe test. In Fig. [7] we show the 
evolution of the piston position (top) and of the gas tem- 
peratures (bottom) with the simulations of the random- 
ized gas previously shown in Fig.[2lJbottom). As one can 
see, the modified (solid lines) and original (dashed lines) 
macroscopic equations generate, a part from ay , two in- 
distinguishable dynamics for x and Ti^r up to times for 
which the mechanical equilibrium is reached. Actually 
for such times the dynamics is essentially given by the 
zeroth order equations (|20ll22p which are the same for 
both the original and modified equations: this explain 
their behavior in the mechanical regime. As one can see 
the main difference in the dynamics is that retaining the 
terms {{V — Vx)v}\ ^ leads the system to stay for a longer 
time interval in a state of approximate mechanical equi- 
librium and thus to a slower relaxation to equilibrium. 
The evolution of the piston temperature May shown in 
Fig. [7] (bottom) clearly show the difference in the two dy- 
namics, and in particular the breaking of equipartition at 
equilibrium for the modified dynamics. 

It should be stressed that the behaviors shown in the 
above figures is not related to the peculiar parameters 
choice as it has been verified in other simulations (not 
shown) . 

The picture emerging from this comparison is that the 
model we introduced goes to the correct equilibrium state 
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FIG. 7: (Color online) (top) Evolution of the piston position 
for parameter as in Fig.[2jbottoni) now plotted in log-lin scale. 
The (red) symbols refer to the randomized gas simulations, 
the dashed curve to the evolution predicted by the macro- 
scopic equations and the solid curve (which perfectly super- 
imposes on the simulation data) is explained in the text, (bot- 
tom) Same as (top) but for the evolution of the gas tempera- 
tures Ti and Tr as in label. The two curved which start from 
zero show the evolution of the piston temperature Tp = May 
as obtained from the macroscopic equations (dashed curve 
which saturates to Tp « 100 and as explained in the text the 



solid one which saturates ad Tr, 



200. 



and respects the basic physics principles as conservation 
of energy and its equipartition. On the other hand, the 
relaxation is slower than the one of the microscopic model 
it should correspond to (i.e. the randomized gas). As 
shown in the above figures the terms responsible for this 
slowing down arc the ones coming from the fluctuation 
of the piston velocity in the expansion of the function G. 
At present wc dot not have a definite understanding ei- 
ther of this difference in the dynamics of the macroscopic 
equations and of the randomized gas, or of the very good 
agreement between the modified equations and the ran- 
domized gas. We suspect that the contribution to the 
relaxation time of these terms might be counterbalanced 
by the resummation of higher order terms in the ^ ex- 



pansion. However, carrying the analysis to higher order 
terms is very demanding due to the proliferation of terms 
in the expansion, and therefore we shall not discuss it 
here. 

We conclude this section by mentioning that the equa- 
tions derived by Gruber and coworkers [l9| predict a rea- 
sonable relaxation time and equipartition at the same 
time. We remark that going to higher order in their ap- 
proach requires to include into the description also higher 
order moments of the piston velocity, namely V^, . . .; 
while within our approach the description remains at the 
level of the second order moment because the statistics 
is constrained to remain Gaussian (Maxwell-Boltzmann) 
for the gases and the piston. This may be the origin of 
the difference between our and their derivation. Their 
equations at the first order might likely correspond to 
a resummation of higher order terms in our approach, 
meaning that their closure with the second moment of 
the velocity is exact while ours remains only approxi- 
mate. 



IV. CONCLUSIONS 

In the framework of kinetic theory, wc derived a set 
of deterministic equations describing the evolution of the 
macroscopic variables in the adiabatic piston problem. 
Our basic assumptions are that at each time the gases 
in the two compartments are perfect, spatially homoge- 
neous, and described by the MaxwcU-Boltzmann statis- 
tics. Thus, at the level of simulations, a (randomized) 
gas model has been introduced with aim to have micro- 
scopic model respecting such assumptions. We obtained 
a set of five ordinary differential equations for the vari- 
ables that describe the macroscopic state of the system, 
namely the mean position of the piston, the average ve- 
locity of the piston, the temperatures of the gases in the 
two compartments and the second moment of the piston 
velocity. The equations are derived up to the first order 
in — 

At the zcroth order they describe a deterministic piston 
characterized by a velocity distribution collapsed on the 
mean, namely P{V) = 6{V — Vx)- This is enough to 
solve the problem of finding the final state of mechanical 
equilibrium and the result coincides with that derived in 
Ref. [131 by using a different approach. 

At the first order the fluctuations of the piston veloc- 
ity, now assumed to be Gaussian, allow for recovering 
the correct final thermodynamic equilibrium. Although 
the evolution of the macroscopic observables provided by 
this set of equations is in good qualitative agreement with 
simulations of the randomized gas, we found some quan- 
titative discrepancy for the relaxation timescales. 

Apart from the performance in comparing with the 
simulations, we would like to stress the conceptual as- 
pects of the method wc developed. It allows for a trans- 
parent description of the macroscopic dynamics of a non- 
trivial non-equilibrium problem, similarly to how the per- 
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feet gas law ean be derived from the microscopie eollisions 
by using elementary kinetic theory. 
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